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Abstract 



> 

QQ I In this paper, we consider a method for implementing a quantum logic gate with 

^^ ' photons whose wave function propagates in a one-dimensional Kerr-nonlinear pho- 

^^ ■ tonic crystal. The photonic crystal causes the incident photons to undergo Bragg re- 

\^ , flection by its periodic structure of dielectric materials and forms the photonic band 

^D I structure, namely, the light dispersion relation. This dispersion relation reduces the 

group velocity of the wave function of the photons, so that it enhances nonlinear 
Mh, interaction of the photons. (Because variation of the group velocity against the 

■y I wave vector is very steep, we have to tune up the wavelength of injected photons 

j^ ■ precisely, however.) If the photonic crystal includes layers of a Kerr medium, we 

^ , . can rotate the phase of the wave function of the incident photons by a large an- 

gle efficiently. We show that we can construct the nonlinear sign-shift (NS) gate 
proposed by Knill, Laflamme, and Milburn (KLM) by this method. Thus, we can 
^ ■ construct the conditional sign-flip gate for two qubits, which is crucial for quantum 

computation. Our NS gate works with probability unity in principle while KLM's 
original one is a nondeterministic gate conditioned on the detection of an auxiliary 
photon. 

1 Introduction 

Since Shor's quantum algorithms for prime factorization and discrete logarithms appeared, 
many researchers have been devoting their efforts to realization of quantum computation 
in their laboratory [1]. To implement a quantum computer, we have to prepare qubits, 
which are two-state systems, and quantum logic gates, which apply unitary transforma- 
tions to qubits. Generally speaking, we can apply an arbitrary transformation of U{2) 
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group to a qubit at ease no matter which physical system we choose as the qubit. How- 
ever, almost all the researchers think implementation of a two-qubit gate to be very 
difficult because the two-qubit gate has to generate entanglement between two local sys- 
tems. Moreover, it is shown that we can construct any unitary transformation applied to 
qubits from U{2) transformations and a certain two-qubit gate, such as the controlled- 
NOT gate or the conditional sign- flip gate [21 [3l [^ . From the above reasons, theoretical 
and experimental physicists aim for realizing the two-qubit gate. 

Among proposals made for implementing quantum computation, for example realiza- 
tion of qubits by cold trapped ions interacting with laser beams [SIE], polarized photons in 
the cavity quantum electrodynamics system [7] , nuclear spins of molecules under nuclear 
magnetic resonance [8], and so on, the scheme of Knill, Lafiamme, and Milburn (KLM) 
is very unique O [10] . They show a method for constructing the conditional sign-ffip gate 
from linear optical elements (single photon sources, beamsplitters, and photodetectors). 
Because many researchers believe that only nonlinear interaction between two qubits can 
generate quantum correlation, KLM's scheme that does not require any nonlinear devices 
seems to be novel and attractive. In KLM's method, Bose-Einstein statistics of photons 
play an important role, so that we can observe bunching of photons in outputs of beam- 
splitters. Recent progress about KLM's linear optical quantum computing is reviewed in 
Ref. dU. 

To implement the conditional sign-flip gate, KLM use the nonlinear sign-shift (NS) 
gate, which applies the following transformation to a superposition of the number states 
of photons |n) for n = 0, 1, 2: a|0) + (3\l) + 7|2) -^ a\0) + (3\1) - 7|2) (with probability 
1/4). The NS gate is essential for KLM's scheme, and they construct this gate only from 
passive linear optics. Although KLM's NS gate is a nondeterministic gate that works with 
probability 1/4, we can detect its false event by measuring an auxiliary photon. 

In this paper, we discuss the method for implementing the NS gate with a one- 
dimensional Kerr-nonlinear photonic crystal. 

A photonic crystal is a periodic structure of dielectric materials whose dielectric con- 
stants are different from each other. If a wavelength of a propagating electromagnetic wave 
is comparable with the period of the crystal, Bragg reflection occurs and the photonic 
band gap (the light dispersion relation) is formed [T2| IT3[ [T3] . This dispersion relation 
reduces the group velocity of the electromagnetic field of incident photons, and thus it 
enhances nonlinear interaction of the photons. If we construct the photonic crystal from 
alternating layers of a Kerr medium and a medium with linear polarization (a medium 
with x^'^^ = and x^^^ 7^ and a medium with x^"^^ = X^^"* = 0? where x''"^ is the nth-order 
nonlinear electric susceptibility) as a one-dimensional periodic structure, we can rotate 
the phase of the wave function of the incident photons by a large angle efficiently. (In 
general, for almost all the dielectric media, Kerr nonlinearity is too weak to rotate the 
phase of wave function of photons by a certain angle.) Using this method, we realize 
KLM's NS gate. In our implementation, the NS gate works with probability unity in 
principle. 

Someone may make an objection to our proposal because we are going to introduce 
a nonlinear device into KLM's scheme that aims to construct a quantum computer only 
with passive linear optics. However, the author thinks that our proposal is a practical 
compromise because the Kerr medium is a well-studied material in the field of optics and 



rapid progress has been made in experimental studies of photonic crystals. The idea of 
making use of the Kerr medium for the conditional sign- flip gate can be found in Ref . [15] . 
Inoue and Aoyagi show experimental evidence of enhancement of the nonlinearity in a 
Kerr-nonlinear photonic crystal in Ref. [IB] . Numerical investigations of a one-dimensional 
Kerr-nonlinear photonic crystal are made in Refs. [TTJ [18]. Some research groups have 
recently proposed optical quantum computation schemes using weak cross-Kerr nonlin- 
earities and strong probe coherent fields [TH [201 [SB [221 [23]. Their proposals and our 
method use optical nonlinearities in common with each other. 

This paper is organized as follows: In the rest of this section, we give a brief review 
of KLM's scheme. In Sec. [21 we outline the NS gate realized by the Kerr-nonlinear 
photonic crystal. And then, we estimate the time-of-flight that the injected photons need 
for obtaining a certain angle of the phase rotation caused by nonlinear interaction in the 
homogeneous Kerr medium. In Sec. [31 we investigate the photonic band structure in the 
photonic crystal and evaluate the group velocity of the photons. From these results, we 
give a concrete example of a design for the photonic crystal that realizes the NS gate 
with giving numerical values of the period of the crystal, the thickness and the dielectric 
constant of each layer, and so on. In Sec. [H we give discussions. In Appendix |Al we 
explain the third-order nonlinear susceptibility x^^^ and the nonlinear refraction coefficient 
n2. In Appendix [Bl we derive the wave equation and the effective Hamiltonian of the 
electromagnetic field in the nonlinear dielectric medium. 

Here, we sketch out KLM's quantum circuit for implementing the conditional sign- 
fiip gate. First, we construct a qubit from a pair of optical paths. The optical path 
(mode) forms a physical system which takes a superposition of the number state \n) for 
n = 0, 1, 2, ..., where n is the number of photons on the path. 10)^.1 (8) |l)a;2 is a state where 
modes xl and x2 have zero and one photons, respectively, and we regard it as a logical 
ket vector |0)a;. We regard |l)i.i ® |0)a;2 as a logical ket vector |i)a;, similarly. And we 
describe an arbitrary state of a qubit as 10)^. = a\0)x + P\i)x for |ap -|- |/3p = 1. (This 
construction of a qubit is called the dual-rail qubit representation [T5].) 

Second, we define the NS operation as follows: 

\^) = a\0)+P\l)+^\2) 

^|V^') = «|0)+/5|1)-7|2). (1) 

We pay attention to the fact that the NS operation does not work on a qubit but on a 
superposition of the number states of a single mode, |0), |1), and |2). 

Third, we construct the conditional sign-fiip gate from the NS operation as shown in 
Fig. m An optical network drawn in Fig. [1] works as the conditional sign-fiip gate, whose 
operation is given by 

for J, A; G {0,1}. (2) 

Let us confirm the function of this network. In Fig. [H symbols Bl and B2 represent 
beamsplitters, which transform the incident number states of modes al and a2 as follows: 
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Figure 1: Implementation of the conditional sign- flip gate with the NS operations. Qubits 
X and y consist of a pair of modes xl and x2 and a pair of modes yl and y2, respectively. 
Symbols Bl and B2 represent beamsplitters. Symbols NSl and A^5'2 represent the NS 
gates. Photons travel from left to right in this network. 
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forn,mG {0,1,2,...}, (3) 

where a\ and al are creation operators of modes al and a2 respectively, and their commu- 



tation relations are given by 
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beamsplitters Bl and B2 replace aj and al with {l/^/2){al + al) and {l/\/2){a\ 
respectively. We give attention to the fact that B2 applies an inverse transformation of 
Bl. Symbols NSl and NS2 represent the NS gates that apply the operation given in 
Eq. ([1]) to the modes al and a2, respectively. 

If we put a superposition of |0)^|0)j,, |0)^|l)j/, and |l)2;|0)j/ into the left side of the 
network shown in Fig. [H the network leaves it untouched and returns it as an output from 
the right side of the network. However, if we put a state |l)2:|l)y = |l)rci|0)x2|l)j/i|0)j;2 
into the network, the following transformation is applied to the modes xl and yl: 
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Thus, we obtain — |l)a;|l)y = — |l)a;i|0)x2|l)j/i|0)j/2 as an output for the input state |l)j;|i)y. 
Therefore, we confirm that the network shown in Fig. [1] realizes the conditional sign- 
fiip gate. Moreover, we notice that implementing the NS gate is essential for the network 
in Fig. [H 



2 An outline of the NS gate 



In this section, we explain the NS gate, which is constructed from a one-dimensional 
Kerr-nonlinear photonic crystal. As shown in Fig. [2], we inject photons whose wave 
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Figure 2: An outline of the NS gate constructed from a one- dimensional Kerr- nonlinear 
photonic crystal. The photonic crystal is composed of alternating layers of materials A 
and B with thicknesses I a and Ib, respectively. [Thus, a period of the crystal is given 
by {I A + ^b)-] The material A is a linear dielectric medium whose dielectric constant is 
given by tA, and the material 5 is a Kerr medium whose dielectric constant and third- 
order nonlinear susceptibility for optical Kerr effect are given by e^ and Xb^ respectively. 
A wave function of injected photons \ip) propagates through the photonic crystal with 
nonlinear phase rotation of the optical Kerr effect enhanced by the slow group velocity, 
and it goes out of the photonic crystal as \ip'). 



function is given by \ip) (a superposition of |0), |1), and |2)) into the photonic crystal. 
The photonic crystal is composed of alternating layers of materials A and B, and they 
form a one-dimensional periodic structure with a period {I a + Ib)- The material A is 
a linear dielectric medium whose dielectric constant is given by eA, and the material B 
is a Kerr medium whose dielectric constant and third-order nonlinear susceptibility for 
optical Kerr effect are given by es and Xb^ respectively. While the wave function of the 
photons travels through the photonic crystal, it undergoes Bragg reflection and its group 
velocity decreases. At the same time, the phase of the wave function is rotated by the 
optical Kerr effect induced by the layers of he Kerr medium B. Because the slow group 
velocity enhances the Kerr nonlinearity, we can expect the phase to rotate by a large angle 
efficiently. 

Here, we estimate the rotational angle of the phase caused by the optical Kerr effect in 
a homogeneous Kerr medium by way of trial. (We do not investigate the phase rotation in 
a one-dimensional Kerr-nonlinear photonic crystal here. We evaluate the phase rotation 
caused by the photonic crystal in Sec. O) According to the Refs. |2H [25], we have the 
effective Hamiltonian of photons in the homogeneous Kerr medium with e and x^^^ as 



n ~ hu{a^a + -) + hxia^fa^, 



and 



X= o" 



9/icj^X 



2^(3) 



(5) 



(6) 



8 t^Sd ' 

where tu is a circular frequency of the photons, S represents the cross-sectional area of a 
device made of the Kerr medium, and d represents the width of the wave-packet of the 
photons in the medium. Thus, Sd is the volume for the quantization of photons. [While 



-^(3) jg given by a fourth-rank tensor as P = x^^^ '■ E + x^^^ '■ EEE in general, x*''^^ iii 
Eq. ([6]) represents a certain component of the fourth-rank tensor. In Appendix \^ we 
give a brief explanation about x^^^ ^^id the nonlinear refraction coefficient 71,2, which is 
observed in an experiment direct. In Appendix [Bl we show derivation of the effective 
Hamiltonian given by Eqs. ([5]) and ([6]).] 

To derive the effective Hamiltonian given by Eqs. ([5]) and ([6]), we carry out the fol- 
lowing process. First of all, we perform the quantization of the free electromagnetic field 
whose Hamiltonian does not include x*^^"*- Next, regarding x*^^) as an expansion parameter, 
we evaluate the first-order perturbed energy, which we can consider to be the effective 
Hamiltonian. (This process is explained in Appendix [Bl) In the derivation of the Hamil- 
tonian given in Eqs. ([5]) and ([6]), we assume that photons travel in the homogeneous 
Kerr medium. However, the system that we discuss in this paper is photons propagating 
through the photonic crystal. Thus, we have to consider the wave function of photons 
which extends over many layers of the photonic crystal. Therefore, the Hamiltonian given 
by Eqs. ([5]) and ([6]) offers us a rough picture of the phase rotation. (Moreover, although 
we discuss the phase rotation by quantum mechanics in this section, we derive the pho- 
tonic band gap by classical theory of electrodynamics in Sec. [31 Hence, our analysis of 
the quantum gate is semi-classical.) 

According to Eq. ([2D, the number state of the photons \n) for n = 0, 1,2, ... obtains 
the phase rotation exp[2xn(n — l)t], where t is photons' time-of-flight in the crystal. 
[We remember (a^)^a^|r2) = n{n — l)\n). Because the first term of Eq. ([5|) (namely, the 
unperturbed Hamiltonian) gives only the background phase rotation exp{za;[n-|- (l/2)]t} 
for In), we can neglect it.] Thus, the Kerr medium transforms the wave function of the 
traveling photons as a|0) + /3|1) + 7|2) -> a|0) + (3\l) + e*2^*7|2). If we let \xt\ be equal 
to 7r/2, we can use this photonic crystal as the NS gate. 

In Refs. [261 [23 [28], nonlinear optical properties of GaAs/GaAlAs multiple quantum 
well (MQW) material are studied. In Ref. [27], Miller et al. prepare a sample that consists 
of 84 periods of 144-A GaAs and 102-A GaAs/Gao.yAlo.aAs and obtain the nonlinear 
refraction coefficient n2 ~ 1.2 x 10^^ [cm^/W] at wavelength (~ 850 nm). The dielectric 
constants of GaAs and GaAlAs are given by e = 13.0 x eo, where eo(= 8.85 x 10~^^ 
[F/m]) is the dielectric constant of vacuum. From the result of ^2 ~ 1.2 x 10~^ [cm^/W], 
the third-order nonlinear susceptibility of GaAs/GaAlAs MQW structure is given by 
X^^^ = 2.97 X 10~2 (esu)= 3.67 x 10"^^ [m-C/V^] (SI unit). This |x^^^| is much larger than 
those of other materials. [The relation between 77-2 and x^'^^ is explained in Appendix [XI 



Let us estimate t = 7r/(2x), the time-of-fiight for realizing the conditional sign-fiip 
operation, numerically. We assume that the wavelength of the incident photons is given 
by A = 8.47 x 10"^ [m] [u = 2tt{c/X) = 2.23 x 10^^ [rad/s], iy = c/X = 3.54 x 10^^ [1/s]] 
which is slightly out of tune from 850 nm. (We use c = 3.00 x 10*^ [m/s].) We also assume 
that the incident photons form a wave-packet whose cross-sectional area and width in 
vacuum are given by S" = 1.0 x 10~^ [m^] (a square with the side-length 0.1 mm) and 
c?o = 4.24 X 10~^ [m] (five times as long as the wavelength), respectively. This wave- 
packet corresponds with a 14.1 femto-second pulse. Because the width of the wave-packet 



is shortened in the medium as d = -i/(eo/io)/(e/^)'^0; where /io and /i represent the magnetic 
permeabilities of vacuum and the material respectively, we obtain d = 1.18 x 10^® [m] in 



GaAs/GaAlAs. (We use the fact yu ~ /ig for GaAs and GaAlAs.) The Planck constant 
is given by h = /i/(27r) = 1.05 x 10~^^ [J's]- Substituting the above numerical values 
into Eq. (^, we obtain x = 1-38 x 10^° [1/s]. Thus, we need t = tt/{2x) = 1.14 x 10~^° 
[s] for the time-of-flight in GaAs/GaAlAs MQW material. To obtain the time-of-flight 
t = 1.14 X 10^^° [s] in the GaAs/GaAlAs MQW material, we have to prepare a device 
whose length is equal to L = ■y/(eo/io)/(e/i)ct = 9.49 x 10~^ [m]. (In this evaluation, 
we assume that the velocity of the photons in the Kerr medium is not affected by the 
term of x- This is because the contribution of the term of x is much smaller than the 
Hamiltonian of the free photons. We notice that the term of x in the Hamiltonian increases 
in proportion to the square of the number of the photons and we consider the case where 
the number of the photons is equal to two at most now. Details of this discussion is given 
in Sec. |3l) 

In the above estimation, the volume for the quantization of photons Sd is crucial. 
We assume S' = 1.0 x 10^^ [m^] and this quantity corresponds with the cross-sectional 
area of the device of the GaAs/GaAlAs MQW material. To prevent the wave function of 
photons from tunneling outside of the box for quantization, we have to wrap up the device 
of GaAs/GaAlAs in a material whose refractive index is larger than refractive indexes of 
both GaAs and GaAlAs. This means that we have to make a cavity of the GaAs/GaAlAs 
MQW material for confining photons inside, and the cross-sectional area and the length 
of the cavity is given by S* = 1.0 x 10^^ [m^] and L = 9.49 x 10^^ [m], respectively. 

The author cannot judge whether this requirement is feasible for realizing in the lab- 
oratory. Because the ratio between two side-lengths of the cavity is very large, it seems 
to be difficult to fabricate this device. To overcome this problem, we reduce the group 
velocity of photons in the Kerr medium by the dispersion relation induced by the photonic 
crystal. 

3 Dispersion relation of the Kerr-nonlinear photonic 
crystal 

In this section, we investigate the photonic band gap (the dispersion relation) induced in 
the photonic crystal of Fig. [2l And then, we derive the group velocity of incident photons 
from this dispersion relation. 

The wave equation of the electromagnetic field in a one- dimensional Kerr-nonlinear 
photonic crystal is given by 

—E{z,t) - fio<z)-^E{z,t) 

~f^oX^'\z)^Eiz,tf = 0, (7) 

l^E{z,t) = -^^B{z,t), (8) 

where 

e{z + Ia + Ib) = e(-2) for — oo < x < cxo, (9) 



. X J eA(> 0) Q<z<Ia 
'^'> \ eB{>0) Ia<z<Ia + Ib ' 


(10) 


X^'^\z + Ia + Ib) = X^'^\z) for -oo < X < oo, 


(11) 


^(3). x_ f Q<Z<Ia 

^ ^'^~\ X^'^ Ia<z<Ia + Ib • 


(12) 



In Eq. (jTj), we assume that the photonic crystal is periodic in the 2;-direction and uniform 
in the x-y plane, and the electromagnetic wave propagates in the positive z-direction. 
Derivation of Eqs. ([7]) and ([H]) is shown in Appendix [Bl 

Here, we assume that the medium B is the GaAs/GaAlAs MQW material and all 
physical quantities are given in Sec. [2l We can estimate the contribution made by the 
term fXox'^^Kz) {d"^ / dt"^) E {z , tf of Eq. ([7D as follows: Because of Eq. (^ and Eqs. (|69l), 
(!7T|) . (!76|) in Appendix [HI we obtain the relation 

X^'^E^-X^'^^ = 4^. (13) 

eSd 9lj 

(We pay attention to the fact that the number of photons in the wave-packet is equal to 
two at most.) And, using x = 1-38 x 10-*^° [1/s] and oo = 2.23 x 10-*^^ [rad/s], we obtain 

-^ = 5.50 X 10-^ (14) 

9uj 

Thus, we arrive at 

X^'^E^ < e. (15) 

Hence, we can neglect the term that include x^^"* in Eq- (0) for evaluating the photonic 
band gap. Moreover, we take monochromatic light as an approximation of the injected 
pulse into the photonic crystal, so that we only need to consider a stationary solution, 

E{z,t) = E{z)e-'''\ (16) 

(In Sec. [21, we assume that the signal injected to the photonic crystal is a femto-second 
pulse of the typical length 5A. Propagation of ultrashort pulses in the Kerr-nonlinear 
photonic crystal is considered in Ref. [ITj. Scalora et al. indicate that the nonlinearity 
causes various effects to the pulse near the band edge. However, because we need hard 
calculations for analyzing the behavior of the ultrashort pulse in the photonic crystal and 
it is beyond the purpose of this paper, we neglect these effects for simplicity. Thus, we 
adopt the approximation by the monochromatic light.) Thus, the wave equation that we 
have to solve is written down as 

where Eqs. ([9]) and (ITOl) are assumed. This problem is known as Kronig-Penney model 
and we can derive its exact solutions l30l [311. 



Let us solve Eqs. Qj, flTO|) . and flTTl) . Solutions of the region I {0 < z < Ia) and the 
region II {Ia < z < Ia + Ib) for Eqs. flTOj) and flT7|) are given by 

Ej{z) = Cjj^^exp{iKjz) + Cj^exp{—iKjz) 

for J =1,11, (18) 



where 
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(19) 



At points where a finite potential step exists, the above solution must be continuous 
together with their derivatives (that is, the continuity conditions). Moreover, because e{z) 
is periodic, we can apply Bloch's theorem to the continuity conditions of the solution. 

Bloch's Theorem: Suppose Ek{z) is a solution of Eqs. (jH]) and (1171) . Then, Ek{z) 
satisfies the following relations. 



Ekiz) = e^'^Ukiz), 

Uk{z + Ia + Ib) = Uk{.z) for — oo < 2; < oo, 

< A; < 



I A + Is 



I A + Is 



(20) 
(21) 
(22) 



Thus, we can require the continuity conditions to Ei{z) and Eii{z) as follows: 

Ei(0) = En(0), 
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dz 
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The above four equations can be rewritten in the form 
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Figure 3: The dispersion relation of the photonic crystal with Ia = Ib = 3.57 x 10 ^ [m], 
^a/^o = 1-0, and e^/eo = 13.0. Broken lines indicates a point where a;(/^+/s)/27rc = 0.843 
(that is, to = 2.23 x 10^^ [rad/s] and A = 8.47 x 10"^ [m]) and {Ia + lB)k = 0.158. This 
point belongs to the fourth conduction band. 



Q = exp[ik{lA + lB)], 
R = exp^iKiils)- 
From Eq. (I27|) . we obtain detM = 0, and it implies 

cos[(/a + lB)k] — cos{IaKi) cos{IbKii) 



(30) 
(31) 



+ 



27^17^11 



sin(/A-f^i) sin(/s-f^ii) = 0. 



(32) 



From Eq. 032p . we can obtain the light dispersion relation uj = uj{k). 

Here, we substitute physical quantities introduced in Sec. [2] into Eq. fl32l) . The im- 
plementation is as follows: The materials A and B are air and GaAs/GaAlAs MQW 
structure, respectively. The dielectric constants of them are given by e^/eo = 1.0 and 
^b/^o = 13.0. We assume the thicknesses of layers that compose the photonic crystal 
are given by Ia = ^b = 3.57 x 10^^ [m]. Using Eqs. fITIJl) and (1521) . we obtain the light 
dispersion relation drawn in Fig. [31 The group velocity of the light in the photonic crystal 
can be written in the form 

^.(fco)=^ • (33) 

fco 

We draw the group velocity Vg{k) of the fourth conduction band in Fig. HI 

Because the variation of Vg against k is very steep in the region where Vg takes a small 
value as shown in Fig. HI we need to tune up the wavelength of injected photons precisely 
for obtaining the slow group velocity. (At the same time, we have to adjust Ia and Ib 
precisely because Vg is sensitive to the thickness of the layers, as well.) 

If we inject photons with A = 8.5 x 10"''' [m] and uj = 2.22 x 10^^ [rad/s] into the 
photonic crystal, its wave vector and group velocity are given by {I a + lB)k = 0.295 and 



10 



Vg/C 




(lA+lB)k 



Figure 4: The group velocity Vg{k)/c of the fourth conduction band in Fig. [3l When we 
have {Ia + Ib)^ = 0.158, we obtain Vg/c = 0.0995. Broken hues indicates this point. 

Vg/c = 0.171, respectively. On the other hand, if we inject photons with A = 8.47 x 10"^ 
[m] and u = 2.23 x 10^^ [rad/s] into the photonic crystal [this corresponds with u!{Ia + 
lB)l'2.'nc = 0.843], its wave vector and group velocity are given by {Ia + Ib)^ = 0.158 and 
Vg/c = 0.0995, respectively. 

Now, we obtain the group velocity Vg of photons in the photonic crystal. However, 
photons injected into the photonic crystal undergo Kerr nonlinear interaction only while 
they travel through the layers of the material B. Thus, we have to estimate the probability 
Pa that the photons stay in the layers of the material A and the probability Pb{= 1 — -P4) 
that the photons stay in the layers of the material B. Because the number of photons in 
a certain box is in proportion to energy contained in it, we obtain the relation 



Pa :P, 



B 



UJ /•27r/w 



/O JIa 

where Ej{z,t) and Bj{z,t) are given by 



/ dtS dz-[eAEi{z,tf + —Bi{z,tf]: 
In Jo JO 2 /io 

^ f2iT/uj I-Ia+Ib 1 o 1 

— / dtS dz -[eBEn{z,tf + —Bn{z,tf] 

Ztt Jo Jia 2 yUo 

""rft r dz[^Ei{z,tf + Bi{z,tf]: 
Jo eoc^ 

/ dt / dz[-^Eii{z,tf + Bii{z,t)% 

Jo JIa enC^ 



(34) 



E,{z,t) 
B,iz,t) 



Re{[Cj+ exp{iKjz) + Cj- exp{—iKjz)] exp(— zcjt)}, 

K 

Re{ — -[Cj-^ exp{iKjz) — Cj^ exp{—iKjz)] exp{—iujt)} 

UJ 

for j = I, II. 



(35) 



[We take an average of the energy over time t in Eq. (IMl) .] 
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Substituting u{Ia + Ib)/'^t^c = 0.843 that corresponds with A = 8.47 x 10"'^ [m] and 
{Ia + ^B)k = 0.158 into the matrix M defined in Eq. f l28|) . we can obtain the coefficients 
of the electromagnetic field, Ci+, Ci_, Cn+, and Cn- form Eq. (1271) . After carrying out 
slightly tough numerical calculations, we obtain 

Pa-.Pb-I: 11.2. (36) 

In Sec. [21 we obtain the time-of-fiight rtof = 7i/{2x) = 1-14 x 10~^° [s] that we need to 
fiip the sign of the phase in the homogeneous GaAs/GaAlAs MQW material. Thus, the 
length of the photonic crystal for realizing the NS gate is given by 

Lp^c = Ttof^^^^^Vg = 3.71 X 10^3 [^]_ (37) 

The photonic crystal with the length Lphc = 3.71 x 10"'^ [m] consists of about 10 400 
layers (about 5200 periods) of materials A and B (the air and the GaAs/GaAlAs MQW 
structure). Moreover, each layer of the GaAs/GaAlAs MQW structure consists of about 
15 periods of GaAs and GaAs/GaAlAs (the width of each quantum well is equal to about 
100-A). 

In the above estimation, we cannot obtain a great reduction in the length of the device 
by the slow group velocity in the photonic crystal. (The lengths of the homogeneous 
sample of the Kerr medium and the sample of the Kerr-nonlinear photonic crystal for 
the phase rotation are equal to 9.49 [mm] and 3.71 [mm], respectively.) To obtain a 
significant reduction, we have to adjust the wavelength of the injected photons more 
precisely. However, as shown in Fig. HI the variation of the group velocity Vg against the 
wave vector k is very steep when Vg takes a small value, so that an accurate adjustment 
of the wavelength of the injected photons is very difficult. This is a weak point of our 
proposition. However, the author thinks that we can expect the great reduction in the 
length of the device by a precise tuning up of the parameters. 

4 Discussion 

According to our estimation, to realize the NS gate, we have to construct the photonic 
crystal from about ten thousands layers, each of whose thickness is equal to about 0.4 /im. 
Moreover, the layers of the Kerr medium has the GaAs/GaAlAs MQW structure. These 
requirements seem to be very severe. However, for example, Noda and his collaborators 
fabricate woodpile-structure of GaAs, whose typical period is equal to 0.7 /xm [321 [33] . 
Hence, the author believes that we can construct the NS gate from the one-dimensional 
Kerr-nonlinear photonic crystal in the not-too-distant future. 

The variation of the group velocity Vg against the wave vector k is very steep in the 
region where Vg takes a small value as shown in Fig. [H Thus, Vg is sensitive to the 
thickness of the layers that compose the photonic crystal and the wavelength of injected 
photons, so that we have to adjust them precisely. In this paper, we show a plan for the 
NS gate with the GaAs/GaAlAs MQW structure as a concrete example. In our example, 
we cannot obtain a significant reduction in the length of the device caused by the slow 
group velocity in the photonic crystal, compared to the homogeneous sample of the Kerr 
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medium. However, the author beheves that careful tuning up of the parameters in the 
laboratory realizes a feasible design of the device. Because the adjustment of the physical 
quantities is very subtle, we may find another design that is better than the plan we have 
shown in this paper. 
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A The third-order nonlinear susceptibihty and the 
nonhnear refraction coefficient 

In this section, we explain the relation between x^^'' (the third-order nonlinear suscepti- 
bility) and n2 (the nonlinear refraction coefficient), which is obtained by an experiment 
for a Kerr medium direct. (Details of this topic can be found in Ref. ||29j.) Using this 
relation, we can calculate x^^^ from an experimental data of n2- 

In general, a nonlinear dielectric polarization of a material is given by a function of 
the electric field as 

j jk jkl 

ioT i,j,k,l e {x,y,z}, (38) 

where x^^'' represents the tensor of the nth-order nonlinear susceptibility. Let us assume 
that an electromagnetic plane wave propagates in the positive z-direction in the material, 
so that the electric field E and the magnetic fiux density B lie in the x-y plane. Letting 
the X axis and the y-axis be parallel to E and B respectively, we can write the electric 
field as E = {E, 0, 0). Moreover, to simplify P defined in Eq. fl38j) . we make the following 
assumptions about x*^"^: 

(2) ^ 

xgL = { ;^^ -L7 '^ - ^ (39) 

and Xij!.. = for n>A. Moreover, we write xi?"* = X^^"*- 

From the above assumptions, we obtain the dielectric polarization P = (P, 0, 0), where 

P = x^^^E + x^^^E^ (40) 

Thus, the electric fiux density is given by D = (D, 0, 0) and 

D = eoE + P={eo + x^'^)E + x^'^E^ (41) 

where eo is the dielectric constant of vacuum. We define the nonlinear dielectric constant 
of the material as 

<E) = ^ = eo + x^'^+X^'^E'- (42) 
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[In general, we regard e(0)[= eo + x^^A as the dielectric constant of the material. 
The refractive index of the matter is given by 



n = /^, (43) 

where /iq and /i are the magnetic permeabilities of vacuum and the material, respectively. 
Assuming |eo + X^^^l ^ |x''^^-E'^| and /i ~ fiQ, which is satisfied for almost all the 
materials, we can expand n of Eq. (H3l) in powers of x^^^ as 

n = no + n^^2 + ..., (44) 

where 



no = 
However, because the nonlinear refraction index n is obtained as a function of the field 




intensity / = (1/2) Je/ fiE'^ in the experiment, we rewrite Eq. (jUj) in the form of 

n = no + 77-2/, (46) 

where 

^(3) rp; ^(3) 

n2 = W-~^^, (47) 

and c = l/^/ioCo is the light velocity in vacuum. 

In Refs. [261 123 [28] , nonlinear optical properties of GaAs/GaAlAs multiple quantum 
well (MQW) material are studied. In Ref. [27j, Miller et al. prepare a sample that consists 
of 84 periods of 144-A GaAs and 102-A GaAs/Gao.rAlo.sAs and obtain the nonlinear 
refraction coefficient n2 ~ 1.2 x 10~^ [cm^/W] at wavelength (~ 850 nm). Substituting 
this result of the experiment into Eq. (H7|) . we obtain x*-^-* = 2.97x 10~^ (esu)= 3.67x 10~^^ 
[m-C/V^] (SI unit). (The value of x^^^ is often represented in cgs-esu unit. We have the 
convenient relation, XesL = 8.1 x 10^^ x Xgi .) In the above calculations, we use the 
following facts: The dielectric constant of vacuum is given by eo = 8.85 x 10^^^ [F/m]. 
Both the dielectric constants of GaAs and GaAlAs are given by e = 13.0 x cq- The 
magnetic permeabilities of GaAs and GaAlAs are given by /i ~ /io, where /xq is the 
magnetic permeability of vacuum. Thus, we obtain uq ~ -\/13.0. The light velocity in 
vacuum is given by c = 3.00 x 10*^ [m/s]. 

The value of Ix^^"*! for GaAs/GaAlAs MQW structure given above is much larger than 
those of other materials. [For example, x^^^ of CS2, which is a typical Kerr medium, is 
given by 3.9 x 10"^^ (esu).] 
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V-D{x,t)=0, 




V ■ B{x,t) = 0, 




V X H{x,t) - -D{x,t) : 


= 0, 


Vx E{x,t) + —B{x,t) = 


= 0. 



B The wave equation and the effective Hamiltonian 
of the electromagnetic field in the nonlinear di- 
electric medium 

In this section, we derive the wave equation and the effective Hamiltonian of the electro- 
magnetic field in the nonlinear dielectric medium. 

First, we consider the wave equation of the electromagnetic field in the nonlinear 
dielectric medium. We use this equation for deriving the photonic band gap in Sec. [31 We 
start from Maxwell's equations in a matter, where there are no currents and no charges, 

(48) 
(49) 

(50) 

(51) 

In these four equations, D and B represent the electric and the magnetic flux densities, 
and E and H represent the electric and the magnetic fields, respectively. Moreover, we 
assume the following relations: 

D{x,t) = eoE{x,t) + P{x,t), 

B{x,t) = fi{x)H{x,t), (52) 

where P represents the dielectric polarization, eo represents the dielectric constant of 
vacuum, and /i represents the magnetic permeability of the material. Here, we assume 
fi{x) = /io, where /iq is the magnetic permeability of vacuum, because almost all the 
materials satisfy it. We also assume that the matter is uniform in the x-y plane. 

Let us consider the electromagnetic wave that propagates in the positive z-direction. 
We can describe E and B as functions of z and t. Because the matter is uniform in the 
x-y plane, we can describe P as a function of z and t, as well. 

Using Eqs. fHSl) . fH9l) . fl50|) . fl5T]) and the above assumptions, we can obtain the follow- 
ing relations: 

|-(eo^. + P.)=0, (53) 

oz 

^B, = 0, (54) 

—By + fxo—{eoE, + P,)=0, (55) 



—B, - fXo—AeoEy + Py) = 0, (56) 



d d 

d-z^^-^'dl^ 

^(eoP. + P.) = 0, (57) 

|^.-|i^.-0, (58) 
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^E^ + |s. = 0. (59) 

|b. = 0. (60) 

From Eqs. fl53|) and fl57|) . we can take 

eo^.(;2,t) + P,(z,t) = 0. (61) 

In a similar way, from Eqs. (I54l) and fl60l) . we can take 

5,(z,t) = 0. (62) 

Here, we assume that the material is a Kerr medium whose nonlinear susceptibilities are 
given by Eq. ( l39l) . Thus, the dielectric polarization is given by 

P, = X^'\z)E, + xf\z)Ef, (63) 

(we assume that x ^-^d Xi do not depend on t), and we obtain 

E,{z,t) = Q. (64) 

Now, E{z,t) and B{z,t) lie in the x-y plane. Letting the x-axis be parallel with 
E{z,t), we can write E{z,t) = {E{z,t), 0,0). Thus, from Eqs. (l56l) and (l58l) . we obtain 
B{z, t) = (0, B{z, t), 0). From Eqs. (155!) and (159!) . we obtain the wave equations for _E(2;, t) 
and B{z,t), 

— E(2;,t)-/ioe(2;)— E(2;,t) 

-/ioX^=^)(z)— E(z,t)3 = 0, (65) 

—E{z,t) = --B{z,t), (66) 

where e(z) = eo + X^^^(-2) and x'''^^ = X^x^- 

Second, we consider the effective Hamiltonian of the electromagnetic field in a ho- 
mogeneous Kerr medium. We assume that the electromagnetic field propagates in the 
positive z-direction, and it is described by Eqs. (IBS]) and fIBB]) . Because the material is 
uniform in all directions, we can regard x^^\ X^^\ and e as constants. 

In general, the energy of the electromagnetic field confined in the box of volume Vol 
is given by 

n= f d^x-{D-E + H-B). (67) 

JVoi 2 

Because we make the following assumptions now, E = {E{z,t), 0,0), D = {eE{z,t) + 

X^^^E{z, t)3, 0,0), B = (0, B{z, t),0), H = (0, {l/fio)B{z, t),0), we can rewrite H given in 

Eq. (jnZD as 

n = Ho + v, 



Ho = I d'x\{eE' + -B% 



V = / d'x^^E^ (6^ 

7voi 2 ^ 
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From now on, we perform the quantization of the free electromagnetic field whose 
Hamiltonian is given by Tio- We regard \^ as a perturbation, and we construct the effective 
Hamiltonian of the quantized field by evaluating the first-order perturbed energy. For the 
quantization of the free field, we write E{z,t) and B{z,t) as products: 

Eiz.t) = -^J2^dt)uiiz), (69) 

B{z,t) = ^^g,(tX(z). (70) 

The free electric field E{z,t) has to satisfy Eq. fl65|) with x^^^ = 0. Assuming the normal- 
ization of Ui{z) in volume Vol = Sd and the boundary condition ^^(0) = ui{d) = 0, we 
obtain 

ui{z) = J— sin— for/ = ±l,±2,.... (71) 

Substituting Eqs. (l69l) . (!70|) . and (ITB into Hq, we obtain 

-Ho = E ^'' (72) 

l=±l,±2,... 

-Hi = liqf + ufqf), (73) 



2 



where 



-, ^ J^^. (74) 

Let us define a new variable pi = qi. Then, we can rewrite Hi as Hi = {l/2){pf + ujfqf), 
and we obtain the relations dHi/dqi = —pi and dHi/dpi = qi. Thus, we can regard qi and 
Pi as canonical conjugate variables. We perform the quantization of the field by taking 
the commutation relation [qi,Pm] = iT^^im- This commutation relation corresponds with 
the following relations: 



qi = Y^("l + "')' C^^) 



Pi = W-^(«J-«z), (76) 

[a/, aj„] = 5im-i and [a^, a™] = [aj , aj„] = 0. Rewriting Hz with the creation and annihilation 
operators a] and ai as 

Hi = hui{a\ai + -), (77) 

we accomplish the quantization of the free field. 

Next, we derive the effective Hamiltonian of the quantized field for V as the first-order 
perturbed energy. We substitute Eqs. (JSID, fTTT]) . and fITBl) into V , and we obtain 

^ = 9^^^^~(^)^ I^ PkPlPmPn 

^ e ^a fc,«,m,n=±l,±2,... 

/■ o kixz Inz mnz rniz 

X / d X sm — — sm — — sm — ;— sm — — . (78) 

Jsd d d d d 
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Using the following formula, 



'^ knz Ittz rrniz nirz 

, dz sm — — sm — — sm — ;— sm — — 
a a a a 

O 

'^k,l+m—n ^k,—l+m—n ^k,l—m—n ' ^k,—l—m—n)j l'"^/ 

we obtain 

^ = j;^^ E PkPlPmPn 

^'^ '^" k,l,m,n=±l,±2,... 

X {—4:Sk,l+m+n + 3Sk+l,m+n + Sk+l+m+n,o) ■ (80) 

Here, we assume that we inject only photons of a certain mode / into the material. 
Thus, terms that contain at {k ^ I) vanish when we apply them to the state vector. 
Moreover, we neglect events where photons of mode k{^ I) are created. Thus, for example, 
we ignore the case where three photons of ui are annihilated and one photon of uj3i{= 3uJi) 
is created. This implies that the effective terms of V cannot contain a|, {k ^ I), as well. 
Hence, the effective potential for photons of mode / is given by 

From some calculations, we obtain 

[a] - ai)^ 
= 6(aJ)^af + Sajai + 1 

+ {aly + 4- 4[{alf + a] + {al)\ + aja^]. (82) 

Here, we concentrate on the terms which conserve the number of photons. Hence, we can 
write the effective Hamiltonian as follows: 

H = huJi{alai + -) 

3 X (^^l\2\at f\2 2 , Q ] ,11 

+ ^^5;^(— ) [6(a]) ai +Sa\ai + l\ 
= Hi + Vi, (83) 



and 



Hi = %uji{a\ai + -) 
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When we estimate the phase rotation induced by Kerr nonhnear interaction, Vi makes 
a main contribution and Tii causes only a background phase rotation. Furthermore, in 
most cases, the relation cui ^ \x^^^\hujf/{e'^Sd) is satisfied and we can take Hi ~ Hi. 
Hence, we obtain the following effective Hamiltonian: 



where 



H = Hi + Vi, (85) 



Hi = huji{a\ai + -), 
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